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SUMMARY 


The equation for steady— state temperature distribution caused by a 
thermal source (or sink) in a flat plate surrounded on either side by 
fluids of different temperature is developed and applications are pre- 
sented in three sections in this report. The applications are as follows: 

(a) A thermal error which occurs when thermocouples are used for 
the measurement of plate temperatures is described, and an 
analytical prediction of this error is obtained. 

(b) A determination is made of the effect on heat— transf er rate 

in pin— fin plates due to the thermal conductivity and thickness 
of the plate metal. 

(c) The temperature-distribution equation is applied to the heat 
meter and a correction factor is obtained which includes con- 
sideration of the effect produced by the flow of heat "around" 
the heat meter. 


INTRODUCTION 


This report describes techniques involved in making certain thermal 
measurements which are necessary in the analysis of heating problems 
in aircraft and contains the general solution for the determination 
of the steady— state temperature distribution caused by thermal sources 
(or sinks) in plates surrounded on either side by fluids of different 
temperature. The terms "source" and "sink" are used to denote a means of 
adding heat to or subtracting it, respectively, from a substance. As 
is shown later, thermocouple leads and fins may be considered to be 
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sources or sinks of heat. The first application of the solution is to 
the calculation of the thermal error of thermocouple temperature 
indications when the thermocouple is employed to measure plate 
temperatures. 

Accurate measurements of surface temperatures are useful in determining 
local heat— transfer rates, for instance, in an exhaust— gas air heater or 
along an airfoil surface and in evaluating thermal stresses and temperature 
distributions which are criterions of the stability and life of a 
heater unit. 

The analytical approach to the determination of a thermal error 
considers the junction of the thermocouple in the plate to be a thermal 
source (or sink depending on whether the thermocouple leads are exposed 
to the hotter or colder fluid). The physical system can be visualized 
best by focusing attention on the case for which the leads are exposed 
to the hotter fluid. The thermocouple leads, in this case, are at a 
higher temperature than the plate because they are not being cooled by 
exposure to a cooler fluid as is the plate. Consequently, heat will 
flow through the leads to the plate thus increasing the temperature of 
the thermocouple junction and the plate in the immediate vicinity of the 
junction, which thus acts as a heat source. The thermal error of the 
thermocouple is defined as the difference in the temperature of the 
thermocouple junction and that of the plate far away (or the equivalent, 
the temperature at any point of the plate in the absence of the 
thermocouple). This thermal error is not to be confused with electrical, 
metallurgical, method of attachment, or other errors to which 
thermocouples are subject. 

The general solution is also applied, in this report, to the 
calculation of the effect of plate thermal conductivity on heat- 
transfer rates for pin— fin plates. Ordinarily the thermal conductivity 
of a thin metallic plate is so large compared with the convective con- 
ductances on either side that the thermal resistance of the plate can 
be neglected in heat— rate calculations. When pin fins are present, 
however, conditions may be obtained that promote heat flow radially in 
the plate from the pin bases. Under these conditions the thermal con- 
ductivity of the plate can become an important factor in some cases. 

The equations which allow calculation of heat rates in pin— fin plates 
of finite resistance are developed and these are compared with the 
usual equations for heat— transfer rates in pin— fin plates which 
postulate infinite thermal conductivity of the plate. 

In a third application the solution is applied to the determina- 
tion of a correction factor for a heat meter when the meter is used to 
measure the heat rate through a plate. The correction factors presented 
before (references 1 and 2) have postulated that all heat flows 
through the plats and the meter in the direction normal to their surfaces. 
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This condition does not result when the thermal conductivity of the plate 
or the plate thickness is large and conditions are such that a part of 
the heat flow occurs parallel to the plate surface into or out of the 
metal directly under the meter. The equation for the correction factor 
presented in this report takes cognizance of this latter condition. 

Experimental verification of the conclusions reached, hy means of the 
analyses presented in this report are lacking. An experimental program is 
being planned, however , which will allow verification of these conclusions. 

This work was conducted at the University of California under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


SYMBOLS 


A 


Au 


b 

c 1? c 2 

f 


L h 


IqU) 


Il(x) 


E 0 (x) 


Ki(x) 

k 


area of heat transfer, square feet 

cross-sectional area of pin fin or thermocouple lead,, square feet 

unfinned area of pin— fin plate, square feet 

distance between in-line pins, feet 

thickness of plate, feet 

constants 

unit thermal convective conductance, Btu/(hr)(sq ft)(°F) 

equivalent unit thermal conductance for flow of fluid over 
insulated thermocouple s, Btu/(hr)(sq ft)(°F) 

unit thermal convective conductance over surface of heat meter,, 

Btu/ (hr) (sq ft)(°F) 

modified Bessel function of the first kind, zero order 

modified Bessel function of the first kind,, first order 

modified Bessel function of the second kind, zero order (for 
a discussion of the rotation used for Bessel functions, 
see appendix) 

modified Bessel function of the second kind, first order 
thermal conductivity, Btu/(hr)(sq ft)(°F/ft) 
equivalent thermal conductivity ^ Btu/(hr)(sq. ft)(°F/ft) 


L 


length of pin fin, feet 
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k 


P 

<10 


r 

r T 

Ec 


T 


t 

t* 


a 

P 

5 


circumferential perimeter of pin fin or thermocouple wire, feet 
heat— transfer rate, Btu/(hr) 

heat— transfer rate through plate in absence of heat meter, 

Btu/ (hr) 

radial distance from source center, feet 
radius of source (or sink), feet 
radius of thermocouple lead, feet 
thermal contact resistance 

temperature of fluid, °F 

temperature of plate, °F 

temperature in pin— fin plate at base of pin, °F 

temperature that plate would obtain in absence of source or at 
infinite distance from source (defined by equation (8)), °F 

equivalent to thermal conductance (defined in equations (14) 
and (15)), Btu/(hr)(°F) 

defined by equation (lO), l/(sq ft) 

thickness of an insulating disk (heat meter) placed over the 

source, feet; also thickness of insulation around thermocouple 
leads, feet 


^defined by q = p~)* (hr)(°F)/Btu 


Subscripts i 

1, 2 two thermocouple leads 

Pl* P 2 two opposite sides of 

s source 

m heat meter 


plate 


ANALYSIS OF FIAT— PLATE TEMPERATURE DISTRIBUTION 


The temperature distribution and heat transfer caused by a thermal 
source or sink on a thin flat plate surrounded on either side by fluids 
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of different temperature can be approximated by an ideal system which is 
defined by the following postulates.* 

(1) That section of the plate occupied by the source (or sink) is 
circular and of infinite thermal conduct ivity (that is, the section is at 
uniform temperature) 

(2) The plate has infinite thermal conductivity in the direction 
normal to its surface (that is, there is no temperature gradient normal 
to the plate surface at any point) 

( 3 ) The temperature of the fluid on each side of the plate is 
uniform and constant (steady state) 

( 4 ) The unit thermal convective conductances of the fluids are 
uniform over the plate surface 

Because it can be shown that the form of the solution will be 
independent of whether the temperature distribution is caused by a source 
or a sink, the distribution will be considered, for convenience, as being 
due to a source. Similarly, one specific side of the plate will be 
considered to be in contact with a hot fluid and the other side in 
contact with a cold fluid, even though the analysis can be carried 
through without knowledge of the direction of the heat flow. 

The solution is obtained by making a heat— rate balance on a 
differential annulus of radius r and width dr which is concentric 
with the source center. (Refer to fig. 1. ) 

A heat balance on the differential annulus consists of the following 
terms : 

The heat flowing in the plate radially from the source into the 
annulus 


q. r = -kA 



— 2jtrbk 



( 1 ) 


The heat flowing in the plate radially from the source leaving the 
outer rim of the annulus 


^r+dr “ ^ (dr) 

r+dr 



( 2 ) 
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The heat flowing into the top of the annulus from the hot fluid 

^1 = f Pl A ( T i - *0 = f Pl 2jtr dr ( T l “ *) (3) 


The heat flowing out of the "bottom of the annulus into the cold fluid 

*3.2 = f p s A (P “ 7 2) = f p 2 2ltr “ T P.) ^ 

The steady— state condition has been postulated so that it is possible 
to equate the heat flowing into the annulus to the heat leaving the 
annulus; thus 


T r + 'll = lr+dr + ^2 


(5) 


Substitution of equations (l) to (4) into equation ( 5 ) gives 


2jtkb 


dr 


( r £) 


dr = 2itr dr 


- t 2 ) - f Pl ( T i - 


Pi' 


( 6 ) 


Rearrangement and simplification of equation (6) gives 


£ Mt 

P 

dr 


+ r 


*p + f 

« - p l P 2 .2 


f T X + f. 


Po 2 


dr 


bk T 


r t 


f + f 
P 1 p 2 


0 


(7) 


The term 


V 1 + V 2 

f p + f p 5 


which will be denoted by too,, represents 


1 2 

the temperature which the plate would attain in the absence of the source,, 

(that is, the temperature of the plate at an infinite distance from the source), 
Thus t — t OT is the rise in temperature of the plate at any point due to 
the presence of the source. 


When a new variable is defined 


V 1 + V 2 

t 1 2 — 

f + f 
P 1 P 2 


= "b — too 


(8) 
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and it is noted that 


dt ? _ dt 
dr dr 


and 


d 2 t f _ d 2 t 
dr 2 dr 2 


the form of equation^ j) can he 


changed giving 


2 

r 


2* s 


d t 


dr c 


+ r 


dt* 

dr 


- 3r 2 t* = 


0 


where 


(3 = 



+ f 


P2 



(9) 


( 10 ) 


This differential equation (9) is a modified Bessel equation for 
which the solution (reference 3) is 



( 11 ) 


In the determination of the constant Cq it can he observed that the 
increase in plate temperature (t* = t — t^), due to the presence of 
the source,, must approach zero as r increases. As a result,, the 
constant Cq must he equal, to zero because the function Iq approaches 
infinity as the argument approaches infinity. Thus the final solution is 



( 12 ) 


or 


t too — 



(13) 


In the evaluation of the constant C2 it is necessary to define 
the source more completely than it is defined "by the conditions given 
in the postulates. The most commonly met source will he described by 
the following conditions: (The simplifying assumptions regarding 

direction of heat flow are retained,) 
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1. Heat flow into the source from the hot fluid can he expressed 
by an equation of the form 

q a = a^ 1 ”! ~ ^g) (ccq_ = Constant) (lM 

where t s is the temperature of the source 

2. The heat flow out of the source into the colder fluid can be 
expressed by an equation of the form 

q.g 2 = a 2 (t s - t 2 ) («2 = Constant) (l 5) 

Equations for the conductances and are presented in 

table I for several physical systems. 


The difference between the heat entering and leaving the surfaces of 
the source disk is the heat conducted into the plate. 


<18! “ *b 2 = -fcA. 




(16) 


This equality becomes 

aqfjq - t 8 ) - 02 (t s - t 2 ) = 2jtr s bkC 2 ^pK 1 (#r a ) (lT) 

Substituting t a = too + t* 8 and solving for C 2 gives 

a lC T l - t J + ^(ip - tj 

Q — A 

K 0 (y$r a )(“i + a^) + 2jrr s bk j/p^ ($r B ') (l8) 


The following material presented will concern itself with' the three 
principal applications of the solutions given here in the order given 
below: 

(a) Determination of thermocouple error 

(b) Heat transfer in pin-fin plates 

(c) Determination of correction factors for heat meters 
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DETERMINATION OP THERMOCOUPLE ERROR 


The determination of plate temperatures by means of thermocouples 
presents a difficult problem whenever the thermocouple leads are 
subjected to a temperature different from that of the plate. Because, 
in general, the leads are at a temperature different from that of the 
plate, it is important that some method be available to estimate the 
resultant error in the temperature indicated by the thermocouple. The 
solution given in this report is easily applied to this purpose for the 
case in which the thermocouple leads are attached in the manner shown 
in f igure (2 ) . 


When it is recalled that 


t— = 
00 


f Pl T l + f p 2 T 2 


Pn 


+ f r 


and the constant C 0 is written out, the solution (equation ( 13 ) ) can 
be written as follows: 


■ts - t m 


ctq( T ^ — too} + 0t 2C T 2 — 


CG*l 4- OCg 



(19) 


The temperature t a is the temperature of the thermocouple 
junction (the source) and thus is the temperature indicated by the 
thermocouple (emf measurement usually by means of a potentiometer), 
whereas too can be considered as the true temperature of the plate. 
The difference in temperature (t s — too} then represents the 
thermocouple error. 


Rewriting the equation slightly and noting that 


JPl 

r J 
p 2 


gives 


'P 


t g - t c 


a l 


1 


— a. 


2 f- 


P„ 


a- 


q + Op + 2 jtr s blri/p 


Kn 


(>K) 


Ko ^ 


( 20 ) 
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Hals equation, although an exact solution of the idealized 
system, is cumbersome and can be simplified with small error by 
introducing a few approximations. 

When the convention is adopted that the thermocouple is on the 
side of the plate for which the variables are denoted by the subscript 1 
the following simplifications can be made in equation (20). In usual 
practice the v.alue of 02 will be much less than a,q (table I, 
systems 4 and 5) and terms including og may be eliminated without 
introducing much error; therefore, 'equation (20 ) can be written 



Two further simplifications can be made. Reference 3 gives asymptotic 
approximations of Kq(x) and E Q (x) as x — >0. 

K-i (x) = — x < 0.05 (22) 

X X 


K q (x) = - log e (|) - 0.577 x < 0.05 (23) 


For the second simplification, ctq can be written as (see table I, 
systems 4 and 5) 


Oq = /f e Pik]Axi + \/f e P 2 k 2 Ax 2 (24) 

For thermocouple leads of equal diameter this conductance reduces to 

Cfq = TtrqVS^ ( + V^) ( 2 5) 

If an equivalent thermal conductivity k Q is defined as 
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tlxe equation can be further simplified, to 


a l = 2rtr l 


( 27 ) 


Insulation 



For the case of a bare thermocouple lead the 
term f e is simply the average unit thermal con- 
ductance f]_ for flow over a small wire. The 
equivalent unit thermal conductance for radiation 
should always be included in fj_. When the leads 
are covered by insulation of thickness 8 and 
thermal conductivity k and when a thermal 
resistance E Q exists between the wire and the 
insulation, the term f a is written to include 
these items as 


Insulated thermocouple 
lead. 

^ + ! + E c A < 28 > 

where A is the heat— transfer area over which the thermal resistance R c 
applies. Combining the foregoing simplifications and noting that for 
equal diameters of the thermocouple leads v&'qi = r g results in the 
following equation: 


Tt - t„ 


1 + 


V 

( fag ^ 

r T /2feE e r T ^log e — 2 0.577/ 


(29) 


For small values of 


with 1 ; therefore 


t — t 

S 00 


T — t 

1 oo 


the term in the denominator is large compared 


T 1 - too 


to ° „ r T 


kpb 


- log e ( ,/|r T j - 0.577 


(30) 
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' t — t 

i g oo 

This equation can he used with small error for values of 


less than 0.10. For larger values the more exact equations (20 ) or (21 ) 
should he used. 


Curves of thermocouple error for four combinations of thermocouple 
lead size and the product of plate thermal conductivity and plate thickness 
are given in figures 3 3 4, and 5. The curves were calculated using 
equation (21 ). 

These curves indicate that the unit thermal convective conductances 
over the plate ^f^ and f^ ^ have small effect on the thermocouple 

error, whereas the error is sensitive to the equivalent unit thermal 
conductance over the thermocouple leads. This fact brings out the great 
importance of proper insulation of the thermocouple leads. It is equally 
apparent that whenever possible, the equivalent thermal conductivity k e 
of the thermocouple leads and the diameter of the leads should be as small 
as possible. 

The following table gives values of k e for several common 
thermocouple leads. These equivalent conductivities are only approximate 
and are calculated for conductivities at room temperature. However, in 
view of the difficulty of accurate determination of the convective 
conductances, and so forth, use of these values of k e at high 
temperatures should cause no hesitation. 


k 


e 


Thermocouple 


Chromel-alumel 
Copper— constantan 
Iron— constantan 
Platinum-point rhodium 


Btu 


(hr)(sq ft)(°F/ft)> 


16 

90 

24 

28 


An estimate of the magnitude of the thermal error encountered when 
an attempt is made to measure the temperature of the surface of an exhaust- 
gas and air heat exchanger may be gained by consideration of the 
following example. 

It is desired to determine the temperature of a 0.030 -inch- thick 
stainless— steel plate exposed on one side to gas at 1500° F flowing at 
150 feet per second and on the other side to air at 100° F flowing at 
100 feet per second. The thermocouple will be attached to the gas side 
of the plate and will be made of Wo. 28 B. & S. gage insulated chromel— 
alumel wire. 
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The unit thermal conductances of the fluids flowing over the plate, 
obtained by using the equations in reference h are 

f p = 16.5 Btu/(hr)(sq ft)(°F) 

f =16.2 Btu/(hr)(sq ft)(°F) 
p 2 

, The unit thermal conductances of the fluids flowing 'over the thermocouple 
leads are determined, however, from the equation, or the graph, in refer- 
ence 5 because, for the preceding case, the value of the Reynolds modulus 
is beyond the range of validity of the equation given in reference 4. 

Use of the graph in reference 5 yields, for the unit thermal conductance 
over the thermocouple leads: 


f 1 =f 2 = 160 Btu/(hr ) (sq ft)(°F) 


The temperature indicated by the thermocouple ordinarily would be 
calculated from equations (20), (21 ), or (30), but the solution of 
equation (21) for the particular system is presented in figure 3. 
Therefore, calculating k 0 fo = 2190 Btu^/fhr) (°F)(ft3) and noting 


that f 


calculating k e f e = 2190 Btu2/(hr)(°F)(ft3) 
+ f - 32.7 Btu/(hr)(sq ft)(°F), the ratio 


too A 


t(oo 


is found to be 0.077 (fig. 3) from which the thermocouple error (t s - O 
is calculated to be 


tg - too = 0 . 077 (t 1 - t^J) 

= 0.077 x 692 = 53 0 F 

The estimated true temperature of the plate is thus 


too = tg - 53° r 

where t g is the temperature indicated by the thermocouple. 

When a thermocouple is attached to a plate and the leads are placed 
in thermal (not electrical) contact with the plate for several inches 
from the junction, the heat that flows in the leads to or from the 
junction, as the case may be, will be very small (due to the reduction 
in temperature gradient along the leads ) and . consequently, the thermal 
error should be expected to be small. However, the disturbance of flow 
over the leads will cause a localized increase of the unit thermal 
convective conductance. This will make the junction hotter or colder 
depending on whether the thermocouple is on the hot or cold side and 
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thus introduce an error of unknown magnitude. This flow disturbance 
may he reduced further hy flattening or embedding the thermocouple wires. 
If the thermocouple leads in contact with the plate are placed parallel 
to and away from the direction of fluid flow, the installation should 
have a smaller thermal error than the case illustrated in figure 2. 

Also, if the convective conductances along the thermocouple leads 
are approximately equal on both sides of the plate (hot and cold sides), 
the error is small providing the leads are attached to opposite sides of 
the plate because the heat added to the surface by the lead on the 
"hot— fluid side" is just balanced by that carried away by the other lead 
attached to the "cold— fluid side" of the plate. These solutions may be 
used for two purposes: 

(1) To estimate the true temperature which would exist if the 
thermocouple were absent from a knowledge of the temperature at 
the point of attachment t s and the other terms in the equation. 
(If a thermocouple is present, the true temperature is given 

by tco, the value far away from the point of attachment ) . 

(2) To predict the error t s — too involved for a proposed experi- 
ment when the surface temperature t M and the other terms are 
estimated. 


Remarks 

In connection with the application of the equation for steady- 
state temperature distribution to the determination of the thermocouple 
error, it may be remarked that in any installation of a thermocouple the 
thermal error can be made small if the following precautions are 
observed: 

(1) The thermocouple leads should be brought out on the side of the 
plate where the unit thermal convective conductance over the leads will 
be a minimum. 

(2) The thermocouple leads should be well insulated thermally down 
to the point of contact with the plate. 

( 3 ) The thermocouple leads should be made of metals having low thermal 
conductivities . 

(h) The thermocouple leads should be of small-diameter wire. 

( 5 ) The thermocouple leads should be embedded in the plate material 
if possible. 

Also, when a thermocouple is employed to measure plate temperatures 
the thermal error will: 

(6) Increase with decreasing thickness and decreasing thermal 
conductivity of the plate metal. 
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( 7 ) 
over the 


Increase with, decreasing unit thermal convective conductances 
plate . 


over 


(8) Increase with increasing unit thermal convective conductance 
the thermocouple leads. 


STEAET— STATE HEAT TRANSFER IN PIN-EIN PLATES 


The question as to the effect of the plate thermal conductivity on the 
heat transferred hy a pin— fin plate (fig. 6) arises frequently. 

Ordinarily , for the case of heat transfer through thin metallic plates, 
the resistance of the plate to heat flow can he safely neglected in 
calculations, hut for pin— fin plates In. which there exists a temperature 
distribution along the surface of the plate due to the finite conductivity 
of the plate, the plate thermal conductivity (or thickness) can he of 
Importance . 


Steady-State Heat— Transfer Equations for Limiting Plate 
Thermal Conductivities 

Heat flow from the hot fluid at temperature t_ into the plate can 
he thought of as following two parallel paths into the plate; one 
through the pin fins into the plate, and the other from the hot fluid 
into the nonpinned portion of the plate. 

In the plate two limiting conditions of plate thermal conductivity 
are possible which will affect the heat flow through the pin— fin plate. 

If the plate thermal conductivity parallel to the surface is postulated 
to he infinite, the temperature of the pin bases and of the nonpinned 
section of the plate must he equal because no temperature gradients could 
exist parallel to the surface. The other limiting condition is 
encountered when the plate thermal conductivity parallel to the surface 
is zero. In this case the pin base temperature, t m in the sketch 
accompanying equation (32), would he independent of and, in general, 
different from the temperature, t n in the same sketch, of the non— 
pinned section of the plate because no heat would he able to flow 
from one section of the plate to another. For both limiting cases 
the heat may he considered to leave the plate hy two parallel paths 
similar to those hy which it entered. 

The accompanying sketch presents the thermal circuit for a pin— fin 
plate of infinite thermal conductivity in the direction parallel to the 
surface. The equation giving the heat rate through this plate 
(equation (31)) can he obtained in the same manner as the equation for 
the current in an analogous electrical circuit can he obtained. 
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Similarly the following sketch gives the thermal circuit for zero- 
plate thermal conductivity parallel to the plate surface, and equation (32) 
gives the heat rate for this condition. 


T 1 



2 


The statement that the thermal conductivity of the plate normal to 
the plate surface is infinite, whereas that parallel to the surface is 
zero, finite, or infinite, is a type of idealization that is often 
necessary in obtaining mathematical solutions. Postulation of infinite 
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thermal conductivity normal to the surface is an acceptable idealization 
because the thermal convective resistances on either side of the plate 
are large compared with the thermal resistance of the plate in the direction 
normal to its surface. It will be seen later that the case of infinite 
plate thermal conductivity usually approximates the heat rate through a 
pin— fin plate closely enough to permit use of equation (31 ) > whereas the 
case of zero thermal conductivity parallel to the surface is an idealiza- 
tion which is never applicable to metallic plates.- 

It is clear that calculation of the heat— transfer rates by equations ( 31 ) 
and (32) will yield two limiting values of heat rate between which the 
heat rate for a plate with finite conductivity must fall. 


Finite Thermal Conductivity of Plate 

In order to determine the effect of finite— plate thermal con- 
ductivity or thickness on heat— transfer rates, an analysis of the heat 
flow through a pin— fin plate with infinite thermal conductivity in 
the direction normal to its surface but with a finite value in the 
direction parallel to the plate surface will be made. 

The analysis can be made in the following manner: Heat flow into 

the pin base. 


Ipin = a l( T l “ O (33) 

Heat flow into the unpinned plate, 

Iplate = 0" 1 ~ ^09) (3^-) 

Heat flow into the plate with n pins , 

1 = f p 1 A u( T i “ "O + na l( T l “ H) - q n (35) 

where q n is the decrease in heat flow from the hot fluid to the exposed 
surface of the plate due to its rise in temperature attributable to the 
presence of the pins* 

The equations defining heat flow and temperature distribution in 
the plate are linear homogeneous differential equations so the temperature 
fields and the heat transfer due to each fin are additive* Thus the 
heat-transfer decrease (q n from equation (35)) is the sum of the 
effects of the individual pins* 
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The effect of a single pin is 

pea 

^-decrease = J 2itrf^(t - too) if 
' r s 


( 36 ) 


"but 


t - t w = CgKo (#r) 

ther ef ore 


_ r a E l (^s) 
'^decrease ~ ^ v ^-p^2 ^ 


( 37 ) 


The effect of n pins is 


Tn = 


2nnr g fp 1 c gE- (|/pr a ) 

a 


Thus the heat transferred by the pin— fin plate is 





( 38 ) 


( 39 ) 


The term C 2 can be determined, as before, by mating a heat 
balance on the pin— fin base. When one pin in the plate is considered, 
the following can be ■written: 

Heat transferred into the pin base by the pin. 


ll = a 1 (, T l “ %) 

Heat transferred from the pin base by the opposite pin (or, for 
an unpinned surface, the heat transferred by the flat plate directly 
under the pin base), 

( 41 ) 


= ” r 2) 
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Heat transferred from the pin "base into the surrounding plate. 


Iplate = (^) = ^2*r a b 


r=r c 




(42) 


or 


Opiate = 2rtr a hhC 2 ^K 1 (i/Pr 8 ) 


( 43 ) 


When a heat balance is made the following equation is obtained: 


c 2 = 


a l( T l ~ *0 + - t h) 

2rtr s bk/j3K 1 (i/Pr g ^) 


(44) 


When Cg 


is substituted into equation (39), there is obtained 


dp 


= f P A( T i _ v) + na i( T i - t b) 


nf. 


Pi 


a l( T l - "tb) + 02(^2 ~ t b) 


(45) 


P 1 P 2 


Determination of the pin base temperature t^ will allow calcula- 
tion of q^. The temperature of the base of a pin fin can be considered 
as the sum of the temperature that the base would have with no pins on 
the plate t** • plus the temperature increase due to the pin covering the 
base and the effects of all other pins. 

For in-line equidistant pins (shown in the following sketch) 
utilization of the fact that the four quadrants surrounding any one 
fin are symmetrical gives 
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t-b “ t oo + C 2 E o(l/P r a) 





Plnr-fin spacing. 


+ K, 


0 


(2i£pa) 



4" e • e 



(46) 


O 

o 

o 

o 

o 

o 


The terms in the series 
represent the temperature increase 
of the pin "base due to the 
temperature increments of all 
other pins. This series is for 
equidistant in-line pins hut 
other series for any symmetrical 
pin arrangements 'can easily he 
obtained. It should he noted 
that the number of pins in an 
actual pin— fin plate has no 
hearing on the number of terms in 
the series hut that the series 
should, in any case, he calculated 
for an infinite number of pins 
(see the following discussion). 
Fortunately the series converges 
rapidly and the number of terms 
given in equation (46) will he 
sufficient for most applications. 


Inserting the value of Cg^ which contains the base temperature t^, 
and rearranging results in 


(tb - too) [2*r a bk|/pK ($r B ) 

a l T l + a 2 r 2 ~ t b( a l + 02 ) 



K, 


: 0 (a/p) 



K 0 (3a|/p) + h 0 ( 4a /p) 


+ K, 


(l^Pa) 



(47) 


which may he solved for t-^. Substitution of t-^ in equation (45) then 
gives the heat rate for a pin— fin plate with a finite value of plate 
conductivity. 
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Discussion 

When conditions exist in a pin— fin plate which promote heat flow in 
the plate radially to or from the pin "bases, the decrease in heat trans- 
ferred due to the finite thermal conductivity of the plate metal may be 
appreciable. Figure 7 shows the effect of plate thermal conductivity kp 
for a pin— fin plate chosen to magnify the effect. 

From this figure it can be seen that if the pin— fin plate were made 
of Inconel (kp = 7 Btu/(hr)(sq ft)(°F/ft)) 0.025— inch thick, the heat 

transferred would be about 20 percent less than that for a pin— fin plate 
of infinite conductivity, whereas a pin— fin plate of copper 0 . 025 — inch 
thick would transfer about 2 percent less than the plate of infinite 
conductivity. 

It can be seen from equations (13) and (l8) that the thermal con- 
ductivity of the plate always appears in the product (bkp). For this 
reason whatever is stated concerning the effect of change of thermal 
conductivity on the heat flow through pin-fin plates also applies to 
the effect of change of plate thickness. 

The mathematical analysis of the pin— fin plate has been written 
for a finite area and number of pins even though at one point it was 
necessary to consider the plate as infinite in extent. The justifica- 
tion for this procedure can be found in the following reasoning. 

If one imagines the accompanying 
figure to be a section of an infinite, 
pin— fin plate, it can be stated that 
no heat flows in the plate across the 
midline AB because the tendency 
for flow in the two directions is 
equal. This is equivalent to stating 
that the plate is cut at this line 
and the cut edges are perfectly 
insulated. The temperature dis- 
tribution and the heat flow in both 
sections would be unaltered by the 
change. Similar reasoning for 
other lines will produce an isolated 
section (with insulated edges ) of 
the infinite pin— fin plate in 
which the temperature and heat transfer 
are unchanged from that of a plate 
of infinite extent. 


A 

O OOOOO 
OOOOOO 
OOOOOO 
OOOOOO 
OOOOOO 

B 
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Remarks 

The following remarks apply to the second application of the 
equation for steady— state temperature distribution: 

1. In most cases it is possible to calculate the heat— transfer rate of 
pin— fin plates by assuming the plate thermal conductivity to be 
infinite . 

2. Thin pin— fin plates of low thermal conductivity should be analyzed 
by the method presented if it is to be established that they do not 
appreciably decrease the heat— transfer rate below that which could 
be obtained by a plate of infinite conductivity. 


DETERMIMTION OF STEADY— STATE CORRECTION FACTORS FOR HEAT METERS 


The heat meter (references 1 and 2) allows the determination of heat 
flow through the meter by measurement of the temperature drop through the 
thermal resistance of the meter. Although the thermal resistance of the 
meter is made as small as possible, it nevertheless adds resistance to 
the thermal circuit to which it is applied, thus altering the heat rate 
that this thermal circuit would have in the absence of the meter. For 
this reason it is necessary to apply a correction factor to the heat 
rate through the meter to obtain the heat rate the thermal circuit 
would have in the absence of the meter. This correction factor is 
defined as the ratio of the heat rate through the plate in the absence 
of the meter to the actual heat rate through the meter (that is, 
measured by the meter). When the notation given in the following 
figure is used, the correction factor is qo/l m * 

Methods of correction have been presented in references 1 and 2. 
These corrections, however, have postulated that the heat flows only in 
the direction normal to the plate and meter surfaces (that is, that 
there exists no heat flow "around" the meter). This postulate is 
closely realized in the meter but, due to the high thermal conductivity 



Heat-flow lines in a heat meter. Heat meter shown with zero 
thermal conductivity parallel to surface so that heat-flow 
lines are parallel therein. 
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of metals, is in error when applied to the metal. plate covered h y the 
heat meter. If , for example, the meter is on the hot side of the 
plate, more heat will flow into the plate than into the meter due to its 
insulating effect. Part of the heat which flowed into the exposed 
plate will then flow in the plate parallel to its surface into the 
section under the meter because of the lower temperature of this section 
produced by the insulating property of the meter. Conversely, if 
the meter is on the cold side of the plate, heat will flow from the secticn 
under the meter into the surrounding plate. It is clear then that a 
correction to the measurement of heat rate by the meter which will include 
the effect produced by the heat that flows through the surface and then 
around the meter is desirable. 

In order to obtain a solution it will be necessary to define the 
system by the following postulates. 

(1) The heat meter is circular. 

(2) The thermal conductivity of the heat meter is zero in the 
direction parallel to its surface. 

( 3 ) Thermal conductivity of the metal plate is infinite in the 
direction normal, to its surface. 

(4 ) Fluid temperatures and unit thermal conductances are uniform 
over the plate and meter surfaces. (Unit thermal conductances 
over the plate and meter may differ. ) 

Two solutions for the correction factor for the meter will be 
obtained. The first solution will require an additional postulate that 
the thermal conductivity of the plate metal under the meter is infinite 
(that is, there are no temperature gradients in the plate under the 
meter). The second solution, which will be more difficult to use but 
will approach the actual system more closely, will be based only on 
postulates (l) to ( 4 ), or, in other words, the temperature distribution 
under the meter will be considered a variable in the second solution. 

The analyses will consider the heat meter to be placed on the side 
of the plate for which the variables are denoted by the subscript 1. 


First Solution 

The heat that flows through a flat plate subjected to fluids of 
different temperature on either side is 



2k 
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If a heat meter is now attached to the plate, the heat that flows 
through the meter is given by equation (lk ) 

%. = a l( T l “ t s) (^9) 


where t s is the temperature of the portion of the metal directly under- 
neath the meter (which is considered uniform for this first solution 
(see following figure)). Also 


When 


a l f n A m k m A m 


+ E c 


( 50 ) 


( T 1 - O = ( T 1 - + (H. - ta) 

and it is observed that the system is the same as the one for which the 
general solution (equation ( 13 )) was obtained, the temperature distribu- 
tion in the metal surrounding the meter can be represented by 

t oo~ t = - € 2 K o(\/P r ) 

and at r = r_, 
s y 

too - t a = -Cr£o(ffir^ 

The value of Cg is determined by means of a heat balance on the 
portion of the plate directly under the meter which has been previously 
obtained in the derivation of equation (l8). The values of aq and ctg 
are given in table I Tinder system 3* 



ib 
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Thus the value of q m is 



When the heat rate through the plate when the meter is absent 
(q Q from equation (48)) is divided by the heat rate through the meter 
(q from equation ( 51 )) the correction factor by which the heat- 
transfer rate indicated by the heat meter must be multiplied is 



recalling that 



This equation can be rewritten as 



or as 
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% 


' 1 + X 
“1 


1 + X 
-p Am 



(5*0 


where X 


2rtr 8 hk^ % (i/fog) 
*0 (i/Pr 8 ) 


and the terms — and 3 — r- 
o- f pAm 


are the 


thermal resistances in specific portions of the thermal circuit. 


Because most heat meters are square, it is necessary to determine 
some equivalent source radius r 8 . This radius can he calculated hy setting 
the meter area equal to the area of an equivalent circle and defining r B 
as (see accompanying sketch) 
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Equation (54) can "be used to show the effects of the two limiting 
values of the plate thermal conductivity parallel to the surface (that 
is, k = 0 and k = »). 

For the use of zero thermal conductivity parallel to the plate 
surface, the value of X is zero, and equation (54) reduces to 



( 56 ) 


Thus, for this case the correction factor is the ratio of the 


thermal resistance 


— + — I , which the heat overcomes when flowing 
a 2 


from the hot fluid through the meter and plate to the cold fluid, to 


that resistance 


1 1 
+ 




f -,V 


which the heat overcomes when flowing 


■*1 - Pp 

from the hot fluid through the plate to the cold fluid without the meter 
in place. 


This result is the same as that given hy equation (10 ) of reference 1 
in which case the flow around the meter was postulated to he zero (that 
is, k in the direction parallel to the surface is zero). A slight 
difference, however, exists between equation ( 5 6 ) and equation (10) of 
reference I in that the thermal resistance of the plate perpendicular to 
the surface is not included in equation ( 56 ). The difference introduced 
hy this omission is negligible for metallic plates hut equation (56) can 
he made to include this resistance hy defining a new term f ? which 
can he calculated from *2 



and can he used in the equations in place of 



(57) 


If the other limiting value of the plate thermal conductivity 
parallel to the surface, k equals infinity, is used, then equation (54) 
becomes 
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_X 1 _ 

$0 _ a l a l 

lm “ ~ l 



( 58 ) 


The correction factor for this case can he interpreted as the ratio 

of the resistance to heat flow — , from the hot fluid thro ugh the meter 

a l ± 

and then into the plate, to the thermal resistance — from the 

V* 

hot fluid directly into the plate without the meter in place. 

Table U presents values of q^/g m obtained from equations (54)* 

( 56 )* and ( 58 ) for zero, finite, and infinite thermal conductivity of 
the plate to which the heat meter is attached. 

Second Solution 

A solution more closely approximating the actual heat flow conditions 
can be obtained by including the temperature distribution parallel to the 
surface in the metal under the meter. Thus if the thermal conductivity 
of the plate is postulated finite in the direction parallel to the plate 
surface, the following equations which give the temperature distribution 
for the plate can be written. 
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The temperature distrihution in the plate under the meter Is 
given by the following equation (see equation (ll)) when 0 < r < r m 


■k “ ' t “m = Cil 0 (l/Pmt*) + CgKo (1/^*) 


(59) 


and, similarly, for the temperature distribution in the plate not 
covered by the meter when r m < r < » 


* - \ = 0 l’ I o( / V-) + C 2%(^v) 


(So) 


•where 


Pin 


jf + f-Q 
m P2 1 


tk p 9 (si ft) 


f + f 

3 = fl ± 1_ 

p bk J (sq. ft) 

P 


fp k m C * Btu/(hr)(sq. ft)(°F) 


A area over which contact resistance R c applies, square feet 
8 heat meter thickness , feet 

too temperature which section of plate under meter would attain if 

thermal conductivity of plate parallel to plate surface were zero, 
or, the equivalent, temperature which plate would attain if meter 
were infinite in extent® Defined as 


b _ Vi + 

b «k f m + f Pg 


F 


too temperature which plate would attain if heat meter were absent; 

p 


Vi + Va 




+ f Y 


> °F 
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At r = 0 it is necessary that the heat flow parallel to the plate 
surface he zero, that is: 


or 




0 


Differentiating equation (59 ), 


£ = (ftr) - C 2 (ftr) (61) 

for r = 0, Sj£ = 0, and ((/^.r) = 0, hut Kq = °°; therefore C 2 

must he zero. Thus equation (59) becomes 

t - tco m = Clio ($>) (62) 

At r = 00 it is necessary that t, the temperature of the plate, 
he equal to t^ ; consequently, (t — t„ ) = 0 at r = <»» Thus at r = oo 

equation (6o) g?ves 


t - = 0 -l‘Io (f- v v) + C 2 «K 0 (^r) = o 


(63) 


At r = I 0 therefore the constant Cq* must he zero. As a 

result, equation (60) becomes 


t - toop = C 2 *K 0 (l^p r ) 


(64.) 



MCA TN No. 1452 


31 


Two remaining boundary conditions can "be stated. First, the 
temperatures given hy equations ( 62 ) and (64) must be equal at r = r m . 

G ll 0 (fern) + ^ = C 2 *K 0 (^r m ) + too p (65) 

Secondly, at r = r m equations ( 62 ) and (64) must give the same heat 
flow parallel to the plate surface or, stated mathematically, the 
derivatives of the equations with respect to length r must be equal. 

Thus at r = r m 


Cl l&l (fe) = -C 2 * ^K 1 (j/^r m ) (66 ) 

The heat that flows through the thermopile section (radius r s ) of 
the heat meter is 


writing 



2 ftrf m (t — dr 


* ~ a * (* - Se) + (H. - h) 


(67) 


(68) 


equation ( 6 7 ) becomes 


Mm 


rtr s 2f m(H ~ T ]_) + 2rtf mCi 


rl 0 ( t/f^r) dr 


= jrr ^f* 

s m 


( too m T l) 


2rtf D1 Cqr s l2_ (jfynFo) 




(69) 


Solving equations (65)and ( 66 ) for the constant C-, and inserting 
its value in equation ( 69 ), there results for q m the equation 




In the absence of the meter, the heat q_ Q that would flow through the same area Jtr g 2 
(area of plate covered by thermopile section of heat meter) is 

i° - ” r = 2f p L ( t - p - T i) 

The correction ratio is thus 



(72) 


S&ll ’°M Ml VOW 
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Noting that 


m 


f + f 

V h fp 


f + f . 


m 


P r 


Pm 


equation (72) can he simplified to 


% 

•3m 



f P-]/ f m 


1 

1 




DM 


0 \vrnrm/ T 

(i^% 

) 


(73) 


Table II presents correction factors calculated by means of 
equations (53)* ( 56 ), ( 58 ), and (73) for the particular system pictured 
in the table. The heat-meter correction factors presented In the first 
three columns are for the case in which the meter Is placed -on one side 
of the plate (the side with high f ), and the remaining three columns 
present the results for the case when the meter is on the opposite side 
of the plate (the side with low • fp). 

It can be observed from the results presented in table II that the 
heat-meter correction factors are considerably larger when the meter is 
placed on the side with the higher fp (that is, lowest thermal 
resistance). The importance of the location of the meter is thus 
apparent. A comparison of the results for the case in which the plate 
thermal conductivity is 100 reveals that the two methods of solution 
yield approximately the same value of qo/qm. This is especially so 
when the correction ratio is less than 1.3. Of course the applicability 
of the heat-meter is questionable when the ratio is larger than 1.3. 

It should be noted that the example chosen here is one that will 
magnify the range of correction factors obtained by the various equations 
and that these factors are larger than those ordinarily encountered. 

It will be noticed that the plate thermal conductivity and the 
plate thickness always appear in the product kpb . Thus if it is said 
that the correction factor increases with increasing plate conductivity, 
it may also be said that it increases with the plate thickness in the 
same proportion. 



MACA TE Ho. 1^52 


3 ^ 


The second solution (equation (73)) shows that the larger the value 
of r m , the smaller will he the correction factor which will approach the 
value given by equation (5 6 ), as a limit, as r m is increased. This 
fact has led to the use of guard rings of the same thickness and 
material around the heat meters. 

In actual application the unit thermal conductance may vary along 
the meter surface because the meter acts as a flat plate in an air 
stream (see figure on page 15 of reference 2). The guard rings 
mentioned are useful in minimizing this variation. 


Remarks 

The following remarks apply to the third application of the 

equation for steady— state temperature distribution: 

1 . Two equations ((53) and (73)) have been derived which are useful in 
the estimation of heat-meter correction factors which include 
consideration of heat flow around the meter. 

2. The use of the second, more exact equation (73) is not warranted 
except for the case of large correction factors (say, qo/ Qm > 1-5 )• 

3 . Table II indicates the importance of the location of the heat meter. 

It is apparent that the meter should be located on the side of the 
plate with the highest thermal resistance (lowest unit thermal 
convective conductance). 

4. Table II also shows the effect of variation of the thermal con- 
ductivity of the plate on this correction factor. For the case in 
which the thermal resistance of the plate in the direction normal 

to its surface is small (that is, metallic plates), these results also 
show the effect of thickness of the plate on which the meter is 
mounted. 

5 . Equations (53) and ( 73 ) indicate the advisability of placing a guard 
ring around the meter. 


Department of Engineering 

University of California 

Berkeley, Calif., August 17, 19^ 



MCA TW Wo. 1452 


35 


APPENDIX 


The notations used for Bessel* s functions by various authors differ 
with the result that there is some confusion in their use. A table of. 
equivalence of symbols for Bessel’s functions is given in reference 3 
(p. 64). The modified Bessel function of the second kind v order 

(K v (x)) used in this report is equal to i?ti v+ ^H v ^^(ix) where Hy^^(ix) 

is the Hankel function given in Jahnke and Emde (reference 6). 


The equations presented in this report have been written where 
possible so that the Bessel functions appear (or can appear) as the 
%(x) 

ratio — t — r~. Prom the foregoing discussion it can be seen that 

K 0 (x) 


,%(x) 

K 0 (x) 


-Hq^^ix) 

iH Q ^(ix) 


or in other words it is unnecessary to convert the Hankel functions to 

& (x) 

Bessel functions when the ratio —s — r is used. 

K 0 (x) 


The following approximations (reference 3) will be useful for small 
values of the argument x. 

Por 0 < x < 0.05 


K q (x) = -log e (|) - 0.5TT 


and when 0 < x < 0.05 


Kl(x) = ~ 

Eq(x) 

K^TxJ 


Figure 8 gives the variation of 


as a function of x. 
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TABLE I 

YAHJES OF a^, a^, AMD r 0 FOR VARIOUS SYSTEMS 
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TABLE n 

COMPARISON OF VALUES OF q 0 /q m 


System A 

System B 

b = 

V = °- s 

r s " °'° 
k m " 0J 

1 « Q fl 

~ fu = 20 6 = f =20 

16 l h ( ^4 Pl 

b = f D = 20 6 = A" 

16 1 P 1 I 64 

? f / 

hr i ^ 

A 1 

1 f = 3 
312 ft P 2 

833 ft 

Btu/(hr ) ( sq f t ) ( °F/f t) 

t 1 t 'f =3. 

f h - 3 1 P2 

r m = 0.212 ft r g = 0.0833 ft' 

k m = 0.1 Btu/(hr)(sq ft) (°F/ft) 

k P 

0 

100 

oo 

0 

100 

OO 

First 

solution 

1.10 

1.50 

1.78 

1.10 

1.11 

1.12 

Second 

solution^ 


1.40 



1.11 



^Values obtained from following equations: 

For zero thermal conductivity, equation (56) 
For finite thermal conductivity, equation (53) 
For infinite thermal conductivity, equation (58) 

2 

Values obtained from equation (73). 
























NACA TN No. 1452 


gl 



Figure 2.- Installation of thermocouple in section of plate. 
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Figure 5.- 


Ve> Bt U 2/(hr)(°F)(ft3) 

Summary of thermocouple error for f p ^ + f p = 20 Btu/(hr)(sq ft)(°F). 
<x 2 « ct.^ • 



Figure 6.- Section of pin -fin plate. Pins on one side only. 





Btu/hr 


DO 



k p , Btu/(hr)(sq ft)(°F/ft) 


Figure 7.- Variation of heat -transfer rate with thermal conductivity 
of a pin-fin plate. Number of pins, 144; free unpinned area, 

0.951 square feet; pin spacing, 1 inch; = 0.14 Btu per hour 
per °F. 



Ki(x) 

Figure 8.- Plot of ratio „ / t as a function of x. 
K 0 (x) 
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